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ABbTHACT 

Several reactor models have hean citod hi literature 
for carrying out the polynierization of caprolact.m into 
nylon 6. Those vary from the simple model of batch reactor 
to the complex system of a scries of continuous stirred tank 
rv."ctors followed finally by a plug flow reactor* In the 
present study, the reo.ctor has been modelled by assuming 
i.;hat the top portion of the reactor behaves li.'ce a stirred 
tank reactor and the rest of the reactor as a plug flow reactor. 
1? ^rt of the water fed to the reactor has been assumed to 
.‘Vaooratc immodicatelj'’. The usuall}'' made assumption of a 
corru?! '.tion between the concentrations of tho linear dimer and 
li e. mnnoiuji' has boon eliminated. The modelling equations 
h iv,,; 1.. ui solved using tl\e moment approach. Irop rtios of tho 
no.l.yw-r product ouch ri.n the voight-avorago molecular weight, 
poly diep'-rsi ty indwx . ,nd rjclt viscosity, which, wore not 
obtained so far, have been evaluated. The validity of the 
model in tho case of the plug flow reactor has been chocked by 
comparing some of tho results with tho experimental data 
th.’it have boon published in the literature. 
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INTEODUCTIOM 


i'Tylon 6, obtained by the polymerization of caprolactam, 
has developed infco a polymeric substance of significant 
commercial importance. Caprolactam can be polymerized in 
three ways; (l) hydrolytic, (2) anionic and (3) cationic. 

Of the three, hydrolytic polymerization is the one that has 
been widely followed for the commercial production of Nylon 6. 

Hydrolytic polymerization consists of polymerizing 
caprolactam in the presence of water. Stabilizers or catalysts 
(usually organic acids) and delustering agents are also 
added. The broad steps of the reaction mechanism are as 
follows: 


1* Rin,ti: Onenin/ry : Caprolactam combine o with water to 
form amino caproic acid, 'v/hich is the' linear monomer. 


HN -( 0 ^ 12 ) 5-0 = 0 + H 2 a ^ H2N-(CH2)5-C00H 
2. Polv-addition ; Caprolactam can then successively add 
onto the acid ei,nd-group, forming linear poljnner chains. 


0 


“ H T 

I -CCH 2 ) 5 - C 3-^OH + H N ^(CH2)^-C =0 


H 


0 

It 


H-4 N -(CHp)p-- c } OH 


3* Poly-condensation ; 'Pt^o. such chains formed by poly- 
addition can then react with each other, resulting in the 
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formation oi an amide links.j^e and a molociile of water. 

H 0 _ 0 

HP- H OH + F [HH -(CH2)5 -’i 5-OH 

0 

^ HrHIf-(CHp)^-C ■5- OH + HpO 
^ m+n 

All the three reactions are reversible. If monobasic 
organic acids are used, they can act as chain length regulators 
also., since they can react with amino end groups. In addition 
to these three main reactions, side reactions leading to the 
formation of cyclic oligomers can also take place. 

COOH C=0 + H^O 

f-- 1 

Of the three main .reactions, ring-opening and poly- 
addition are the reactions by which conversion of caprolactam 
into the polyuicr tcHces place and the poly condensation is the 
reaction which determines the final degree of polymerization. 

Once the desired degree of polymerization is achieved, 
any monomer present in the polymer mixture, along with cyclic 
oligomers, is removed either by hot water extraction or by 
applying vacuum. One of the important factors in the process 
is the amount of water fed to the reactor. Large amounts of 
water remaining in the system aid in the reversible reaction 
of the polymer molecule splitting up into smaller units through 
the condensation reactipn. Thus water must be present in 
large quantities initially, but must be removed from the 
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system, once a reasonable degree of conversion of lactam 
has taken place. The optimal conversion at which water 
must be removed can therefore be determined only when a good 
model of the reactor is available. The focus of the present 
work is to develop a mathematical model of the polymerization 
reactors. 

Both continuous and batch reactors are used industrially. 
In the case of the continuous reactor, caprolactam and water 
are admitted at the top of the reactor continuously. The 
upper part of the reactor is normally at about 200°C. The 
reaction mass passes then into the lower region of the reactor, 
which is heated to 260°C. During the passage through this 
zone, water in the reactor mixture starts boiling and some 
of Uio water is removed in the vapour form through the top 
of the reactor. The reaction mass then enters a second 
reactor (consisting of two concentric tubes) at the bottom 
of the ou tier tube and rises up this tube. During this upward 
movement, both residual water and that formed by poly conden- 
sation can escape more easily. The rising water vapour bubbles 
generally cause fixing in the reactor. At the top, the polymer 
melt flows into the inner tube, thereby forming a thin liquor 
film from which any remaining water easily vaporizes. In 
this part of the reactor, boiling does not occur. 

In the case of the batch reactor, the reaction is carro-^d 
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noraally in two stagoo. In the first stage, during which 
ri opcjiing and poly addition reactions occur, the system is 
kept und'.r olovatod press^ires (ohout 4 atms) . After the 
desired degree of conversion of lactam, has been reached, 
vacuum is applied to the reactor, thus removing water. After 
a definite time, vacuum is broken and the polymer melt is 
extruded under pressure. 
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LITiji-aATUllS SURVEY 

Polymerization of caprQlactam , for the commercial 
production of nylon 6, was first started in Geniany in 1940. 

The initial attempts to understand the kinetics of the 
polymerization were made simultaneously by Hermans et al.fll 
.and Kruissink et al. [2] in Netherlands and by Wiloth [3] 
in Germany. While Kruissink et al.. [2] studied the poly- 
merization initiated by tamino caproic acid, Hermans et.al [ll 
investigated the ■pol3rmeriza fcion in the presence of water. 

They concluded beyond doubt that the polynerization process 
occurod through the three b.^sic steps (described earlier), 
a view which was supported by Wiloth [ 3 ]^ 

Wi.'c. •■ai'. ob al fl] contacted batch experiments at 221.5^0 
for bhroe different initial compositions of water and 
caprolactam in the feed, and at 254.5*^0 for one of the three 
compositions. Kruissink et al. [21 who did similar experiments, 
and Hermans et al, "[1] concluded that all the three reactions 
wero catalyzed by the acid end groups of polsaaer chains. The 
rate constants, according to these groups of workers therefore 
consisted of two parts, one being for the uncatalyzed reaction 
and the other for the , reactions catalyzed by the acid end 
groups. The general form of the rate constants was given as 

^The article by Wiloth is in German, hence his I'esults which 
are reported .here are from English sources whei-e this article 
has been referred to. 
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k. = k! + ]c.*' C, where -fclie k.’ and k!’ refer to the nncatalyzod 
and catalyzed parts respectively and C refers to the coheen- 
tration of the acid end groups. Numerical values for the 
kfs and k^s as well as for the equilibrium constants at 
221. 5°C were given by Hermans et al. [1]. In their work, 

Hermans ot al, [1] assirated that the formation of cyclic 
oligomers was low and that their effect on tho kinetics was 
negligible. They also assumed that the concentration of 
the linear dimer was equr.l to that of the linoar monomer. 

kinetic data have also been published by Hoftyzer et al.T^]. 
They used a similar general form for the reaction rate constants 
and gave numerical values for tho rate constants as well as 
equilibrium constants, for a temperature of 260®C. 

fioimschuossol e b ■■'1. fSl carried out the kinetic studies 
in tho ruiigo of 220-265°C and dobormined the temperature 
dcpendonco of the rate constants. The temperature dependence 
was correlated by tho usual Arrhenius relationship. They 
have given the values for the activation energy and the pre- 
exponential factor, for each of the rate constants. 

Among the three sets of rate constants tho.t have been 
given in the literature, the ones estimated by Hermans et al.[l] 
and those by Reimschuessel ot al. f 6] conform with each other, 
so far as the ring-opening and poly addition reactions are 
concerned. Reimschuessel et al. [6] have expressed the 
equilibrium constahts in terms of the heat of reaction and 
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entropy chn,n;^o. i^hcso v.iluos approximate well with those 
■yivoii hy ■'’.er-'v.-.e.':; et al. [ll. Tho various sets ol rate 
constants and equilibrium constants as obtained by the 
different groups of workers have been given in Table I. Here, 
it is pertinent to mention the work of &erdes ot al. [7]. 

Gerdos ct al. [7] determined the rate constants for the 
thre-e basic reactions from tho experimental results of 
Hermans ot al, [1] by a parameter estimation technique. The 
nujinerica.1 values for tho rate constants have not been given. 
However, they appear to claim that tho values given hy Hermiuis 
ot al.[l] are incorrect and that based on the correct values, 
which were not given, considerable improvement in an existing 
process was noticed. In describing the kinetics of the 
process, they have assumed that the concentrations of the 
linear dimer and the linear monomer are related according to 
the chain length distribution of Flory [4]. Howovor, it is 
to bo noted that Flory' s distribution is valid only for 
irrovcrsiblo condensation polymerization. 

All the workers mentioned so far neglected the effect 
of cylcic oligomer formation on the kinetics of the system. 
Mochizuki and 'ito [8] studied tho kinetics of the termination 
reaction (i.e.'the reaction of an organic acid, usually added 
as a stabilizer, with the amino end groups of tho polymer 
chain) and also the formation of the cyclic oligomers. Data, 

i; 

and equilibrium constants and activa.tion energies for these 
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two roo.ctions havo be-an ftivon fSl. The oligomer content at 
vorioua times during the course of polyuviriaation was studied. 
According to their results [8], the concentration ,, of the 
overa,ll cyclic oligomers is never greater than 2 per cent, 

Uork on the modelling of nylon 6 reactors was first 
stirtod by Hoftyzer ot al. [5]- Their aim was to optimize 
the course of temperature and water content in such o, way 
that a curtain DP (degree of polymerization) was obtained 
in the minimum possible time. To simplify matters, they 
presumed tha.t there wore no operational limitations for 
roilizing the existence of such an optimil course of temperature 
liid water content for the pol 3 rmerization process. They also 
acriuiijod batch oporation (or ideal plug flow). By employing 
tiyn Jilt, pro, , 3 ’ ’im; i I'i. ■ , they concluded that the process should 
b»i c. rriud mil; in two stiigosi the first one operating at a 
high content to facilitate the conversion of caprolactam 

and fchu second one at a low water content to attain the 
maximum DP through the condensation reaction. In between the 
two stages water should be removed as rapidly as possible. 

•n-ro limiting values were assumed for the temperature; the 
upper limit being the boiling point of lactam at atmsopheric 
pressure and the lower one corresponding to the melting point 
of the polymer. A degree of polymerization of' 140 was reached 
by this method in the minimal possible time. On- of thu 
drawbacks of this model is that it is rather difficult to 
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iCMlisu Iho i .IS 0 s foiiiov '■1 of prosciit in the 

oysbt-3D li Cur iliu first stn-go, though thoru exo procossos 
uhiCii opu.r.'^/cj on the b-'sic principles upon whichthis 'study 

I’Ms biGuci. fot nhother drawbjick; is that distribution gets 
iiarro\.’cd down. 

'’-i 3 ''‘-subninanian oaid dciiaschuossol f9] analysed the 
aiodellini'.; of ruictors foi' the pol 3 aiierizn,tion process. They 
devised r'. system of a series of perfectly mixed stirred taiilcs 


iollovjed by a plug flow reactoz'. In order to obtain the 
highest possible rates of conversion, they assumed th'-t the 
reactor system consisting of the stirred tanks 1 to (N-l) was 
UGSenti'illy a closed system . Romoval of water was done only 
in the i'Tth reactor, and th.''. t boo at a f ini be rr.', bo instead of 
bh»’ iTiatan baa->oiis removal as w is done by Hoftyzor ot al. [51. 

'.’hi ;.} A' llt.- W'.e do ivrrm.i ib-'d frotu o.n oqo vbion convioc bing the water 
coruu^n Lr . lion ind ch'iin l..‘n;,th concon bration in the Nth 
r.jrictor. I’hcy [9] anal /sod a system consisting o.f throe 
stirred tank reactors in series, followed by a plug; flow 
re/'.ctor, fh'cy found th.at in such a system, considerable 
reduction in the ovur.all reaction time could be obtained, over 
blio plug flow reactor in whj.ch water is removed instantaneously, 
for atbainin.g the sene degree of polymoriz.ation. Such a 
reactor model would bo officienb only when threj or more stirred 
tank reactors in series arc considered. ¥hon only one or 
two stirred tanks are used in series, the sibuation would be 
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difiOront. MoroovoT, they f 9] also have ass'Jia^d that the 
conC'-ntrations of tho linear dimer and linear monomer are 
oqu'il. 

Jacobs o,nd Sch¥eigman [10] have made a detailed study 
of the modelling of the reactors for the pol^nnerization 
process. In their model of the reactor, they have considered 
the reactor to consist of two parts, the first one approximating 
stirred tank reactor due to the strong mixing effect caused 
by bubbling of water vapor and the second part behaving like 
a plug flov; reactor, which is also reasonable because of the 
high viscosity of the polymer at this stage and the presence 
of baffles in this part of the reactor. While formulating 
their ma-thcmatical model, they have based their expression for 
hho lineir dimer on rior;r»o distribution, i.o. , the same that 
l,L’Ui bc„)n by Gordva* ot .al. f'^]. Jacobs and Schweigman 

[10] deteimiincd the oL'fect oi water on the DP at various 
bijmpujL atui'os. From their results it is quite evident that 
uvon a slight variation in the water content in the reaction 
mixture can have a groat influence on the DP. They measured 
the water content in the reaction mixture based on vapour flow 
rate measurements and viscosity measurements for various 
temperatures. On tho basis of the results thus obtained, 
they g.ive an empirical relationship for the determination of 
i^oight fraction of water in the molt at the top of the reactor, 
as a function of temperature. They considered the extent to 
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which \-;ator should be r^jiaovod before poly co:.ic*.uns.ntion 
reaction took effect . They concluded that the increase in 
the relative viscosity wrs hi.^her the greato-r being the per 
c.- nc reraovnl of WvUter from the system. Polymerization at 
elevated pressures with intermediate wntor removal has also 
bec.n considered by them. They found that a 66 oer cent 
decrease in the reaction time could be rachiovod using this 
method over the conventional 7K (vertical tube) column 
I’e 'Ctor, but at the expense of simplicity in the reactor 
decign. 

Considering the various reactor models so far seen, 
i 1 . is evident a reactor based on plug flow assumption alone 
would not be nn adequate one. Also, having seve-.’al stirred 
(i'.iik c , ’.ctoro in series without water removal is unrealistic, 
!.liou ji .hiy conditions dictate a combination of stirred 
tank rutcloKs and plug flct'. reactor. Hence for simplicity, 
the model envifj''.god by Jacobs and Schweigman [10] seems 
plausible. This model has been adopted in the present study. 

The present study aims at calculating properties that 
have not been calculated hitherto. Of all the workers 
referred to, Jacobs and Schweigman [10] 'have calculated the 
viscosity of a 1 per cent solution of nylon 6 in formic acid. 
In the present study the melt viscosity, based on the weight- 
average molecular weight is calculated. Also, all the workers 



12 


'.0 :i-.r h- 
diaur frii 
by B'lory 


VO '■•,cjB’..racd th;'t the concontrations o 
i linear monomer arc equal or relo.tod 
's distribution. The present work ha 
this assumption. 


f the linear 
to e,ach other 
s attempted to 



O I 

o . 


15 


1 


^ I !>• 

O i 
+ ! H 

in 

!>■ i 


o 

C\J 

o 

! H 

o 

CO 

•f 


*«sf- 

m 

m 

. 

1 


+ s ^ 

8 ^ 
i m 


:) I 

HP 

c,i o i a o 

fnO ,^0 

IT\ I'H Ln 


, u. > 

''H ♦ 

U • 

O O r~} 

rH 

'H CM 

•H CM 

’P CM 

!H CM 

O 

*r"| 

'r' •’» t t -P 


fP 


; cri w 

'•Iv- 


<]> 1 ■ 

>-v 

rJ O J a O 


O 

o o : 

+=> o 

dvo I a 


^ • 

u . 

So 1 

1 Uo 

S -P 

t 1 I 

Uo 

OH : 

I f— 1 

) 0 0 , 

, rQ 0 

0 cj ! 

[ rQ 0 

H CM , 

, *H CM 

•H^ i 

1 -HMD 

•H ^ 

•H KO 

-P CM i 

1 H CM 

1 +=> C\J 

H CM 

«P 1 

1 H CM 

O 

•H 

G 

, -H 

0 

•H 

CJ "P 1 

1 ^ -p 

3 s -p 1 

1 -P 

1 cr] 1 

1 P-P 

O c3 , 

O’ rj 


0 ^ cl 

i 0 

j-j 

p.1 ^ I 

t'PqN_. : 

1 Pi'--' 1 

I Pq w 

! ! 

! 


1 


1 


I 



CHAPTJilR 3 


F0RI>!IILATI0N of MODBLLIIG EQUATION'S 


As staged in Chaptor 1, during the pol 3 niierization of 
caprolactam to nylon 6 three basic reactions have been 
r 0 cognized . They are : 


(1) Ring Opening: 


M + W 


S. 


It ! 
^1 


(a) 


(2) Poly Cond onsation: . 


S„ + ^ 

11 m j 




S„_^ + W (t) 

n+ffl ' 

m,n = 1,2, . . 


t . 


3) t’oly Addition: 


(3,1) 


+ M r 

n ^ 




"n+1 

3(X'*~iX- Ip 2 tt • fc 


(c) 


for the symbols used, rofefenco is made to the list 
of s^Tabols, 

Equations for g Plug Flovr Reactor or a Batch Reactor : 

Prom the above three reactions, rate equations for 
the ccfuccntration of ■ caprolactam, acid .and groups, amide 
groups and amino caprotc acid are given fl] as 
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dii 

.'JIT 

dt 

d?. 


- k^(i'iw - IC^S^) - k^(ac - IC^CC-S^)) (a) 




(3.2) 


jf = k ^ Cc '- KgWZ ) + k ^( MC - K3 ( C - S ^)) 


( b ) 

(c) 


dS . 


lihv 


h ‘7 


i . “ i , , 

o 


g~^— = k^(HU-ICj_S^) - 2k2(CSj_-K2W(C-S^)-) 

~k^(MS^-K3S2) (d) 

Tri tho last' equation (i.e. 3.2 d) , ’S 2 ’ rofers to 
conc’-n oration of tho linear dimer. This equation explains 
to vihy a relationship bwtwoon S 2 and had to he assumed 
i,li > previous worJoerc. To avoid this difficulty and to 
•' til Wv; i.idit-avu-ro.go molecular weight the concentrations 
.1.1. .1^.^ n’,...'rled. The ratu expr<..ssions lor tho concen- 

. i-juM.- ot cm bo written .as 


d.l 

rr 


- k ^ 




n=l 


at 


k ^( MW - K3 _ S ^) - 2k2<S 



m=l 

. - K5S2 ) 



(a) 

11=1 


€0 

(b) 


m=l 

(3.3) 


dS , 

n 

rt *”™ 




m 


n-1 

E3 

_, m=n+l m=l 

■kjK^W S^C)a-l)+lC3M(S^^^-Sj^)+kjK3(Sjj^3_-Sj^) (c) 

n 2 
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dorivinc tli.. last aquation in (3.3), it has been 
ai7i:a/.v.-a tha I', all tho (n-l) amido bonds that aro present in 
chainof n units, ax-o equally accossiblo, for the 
x'\,v r -'action. 

'C’ thu' concentration of the acid end groups is 
nothing but tlij sum total of tho concentrations of all 
,.h-. lineal’ ch'iins. Hence 'C* can be denoted by 

CO 

C = ^ (5.4) 

n=l 



To obtain the expression on the right-h;and side of 
‘7qn.(3.7), e '.ch term on the right-hand side of oqn. (3.4(c)) 
ha.^ to be multiplied by (n-l) and the resulting oxprescions 
th-.n aUifliued up. When this is done, oqn. (3. 2 (b) ) is obtained. 
Th.-o". results prove tho consistency of tho present model with 
til', previous rate expressions. It might he poiiiuod out that. 
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i:i t‘.. ^ tion oi tho equation for requires 
rh'ic intv^7r'ition of succjssivn; equations, which is possible 
for irr*;v... rsiblw pol5nnv, I’izations, can not bo carried out 


Tii--.' thruw equations of (3.5) are now made use of to 
;:he Doir.eni: equations. The need for the raoment 
cqti”, I'ioiia will bu givuii in tho next chapter. 

In the kth moment is defined as^ 




i: 

n=l 


c- 

n b, 


n 


(5.8) 


To f’let the kth moment, the expression on the right- 
' : uiue of equation (3.3(c)) is multiplied by n'^ and then 
..ru --iv..''. up from n=2 to n= This is then added to eqn. (3.3(b) ) * 

fh V ar.:! I i r. 


u X 




\ ) -Skg ( ( C-S3_ ) ) -k^ (MS^-k^S2) 

“ - ( Ak4.1-2'"'^S2-M 


+ k 


2 


r 


kl 


T Ax Xk-x + 


XaO 

(Ak-sS-sp 


rn^xTT Ax Ak^x ^2^2* 

k-1 k 


2 2 




^ k4-l~" ^2_ 

k+1 


,k-2t+l \ 
k-2t+l“^ ^2“^1 J 


t=l 


- 2 C\-z'' 5 ^,S,) 

/_-2(C-S^-S2) -(Xk+l-APl -^^2 O(Xt-Si) 


+k3M [j- jirtaji' ^k-j “ (Ak-S3_^^ 

j=sO 



k 


k. 

+v V r T — (-1)^ kj \ \ (-l)kJ 

+.c3.^r ^ tiTk^TT' A k-j - ^ jiTF-TTi 

3=0 . i=0 .. ■- . • ■ 


■'I 




(3.9) 


3=0 


As should bo expected, from the 2nd momonb onwo.rds, !.■ 

X, 2 onwards, the calculation of each moment requires a value 

of the next higher moment. 

The above system of equations, either in 'bhe moment 

dS 

form or direct form (i.e, jtjX) can be used to predict the 
porf orraanco of a plug flow reactor. This part is discussed 
at length in the next chapter. 


liquation s for a Continuous Stirred Tank Reactor ; 

As disceissed earlier-, the CBTR is relevant to model 
bhe top part of the roa.ctor. Both caprolactaiii and water are 
fed to the reactor at roca temperature and this mixture then 
boils duo to the high 'beraperat^ir . at some point in the reactor. 
The reactor upto this stage then is assumed to behave as a 
continuous stirred tank- reactor . Part of the water fed, 
evaporates out of the reactor, while any caprolactjun that 
evaporates, is recycled bade to the reactor, so that the 
influence of the evaporation of caprolactam on the kinobics 
can be neglected. ‘ Clearly this is an idealization of the 
complicated two phase reactor. 

Jacobs and Schweigman [10] have given an empirical 
correlation between the per cent by weight of water in the 


19 


ronction mixture, at tlia top of the molt, and tho 
tempo raturo of tho reactor in this zone. This correlation 
is 

¥ = 1.76 - 0.0060 Tq (3.10) 

where T^ is tho tompero-turo of the melt in this zone. Since 
perfect mixing is being assumed, ¥ then is also tho water 
content in the exit streamof tho CSTfi. 

Whon tho evaporabod caprolactam is rGcj'-cJ.cd, the 
polymerization process in tho CSTR in which tho tomporo.ture 
is loiown, con be describud bjr the following equations; 

- F = 0 (1) 

P.H. - PM + V r^ = 0 (2) 

1 1 M 

(3.11) 

P.C, - PC + rn = 0 (3) 

The subscript ’i' in the above equations refers to 
the inlet conditions, i'-y is the flow of water va-pour leaving 
the reactor, and P is the output flow rate of the reaction 
mass. The 'r’s refer to the rate equations for the respective 
snecies and are given by the right hand sides of equations 

p Tr 

(3.2) and (3.3). Per a CSTR, since = r , the residence 

time, dividing the equations (3.11.2) to (3.11.4) by fY, 
we get 
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I'i 

X 

(J 


i ^ 

T ^ 


-i _ 5 + V. 

T T 


0 


0 


Hi 


n 


r = 0 


( 1 ) 

( 2 ) 

(3) 


( 3 . 12 ) 


% X 

SiHCO th 6 iiilot concontration of poljnnor chain units 
is zero the last equation, i.e. ( 5 . 12 . 3 ) becomes 


n 

X 


r S 


= k 


'n 

n -1 

•2 5 ~ 


<50 


</-> 


Sv, ™ 
n-m m 


+ 2k K.¥ y~ S„- 2 k„ S„ \ S„ 
2 2 m 2 n m 

m=n+l m=]L 


-kjKgW S^(n-l) + k: 3 MCSn_i-S^)+kjK 3 (S^^^-S^) 


n >2 


(3.13) 

k 


i*ultiplyin<;, both sides of the above equation by n 
and DiajTwifiriGe from n=l to :0 yields on the left-hand 

uide and the right hand side will be' the same as the right 
hand side of eqn.(3.9). Prom this, the individual moment 
equations may be obtained by substituting the corresponding 
values for 'k*. For exanple; 

-Aa k^(MU - K^sp - - K2¥(X^- Ap) 

- ( 1 ) 

Xi X \ (3.14) 

—A == k^(]y[¥ - K^S^) + ^ 3 (^X 0 - 

( 2 ) 


« 
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and sixUilar e::pressions may be wri'iten for other moments 
also. 

Substituting the expression for r^^^ the corresponding 
r.'.te equation, eqn. (3.12.1) becomes; 

I = ^ Ao-Si)) 

(3.15) 

The system of non-linear algebraic equations (3.14) 
and (3.15) have to be solved by a suitable technique to 
yield the. molecular weight distribution for the case of CSTR. 

It is to bo noted that still has to be obtained and 
c:iicul.ation of Xj^. still requires Xk+1* 

CSTR I’ollo w od by a Plug Flow Reactor ; 

Thin is nothing but a combination of the mathematical 
jAodeln writ;;on individually for the two cases in the preceding 
pa/'ho, for this combination, the equations for the CSTR are 
solved first. Values of the concentrations of caprolactam, acid 
end groups and also the moments are then taken as the initaal 
values for the plug-flow reactor. The corresponding equations 
for the plug flow reactor are then solved using these initial 
values and a complete picture of the distribution and concen- 
trations of the various species are obtained. 


,,L I. 

TZliO 



ir Vi’iiOD Oi? noLu -, eiogF of ti-ie modelling 
bquatiohs 

In the previous chapter, the formulation of modelling 
equations, in terms of moments, for the p-olymerization 
process, in a CSTR followed hy a plug flow reactor was done. 

It was pointed out that the solution of these equations 
needi.'d a method for the evaluation of the advanced moments 
and the concentration of tho amino caproic acid,S^. 

It io possible to integrate tho rate equations for 

concGiitra cions of the chain lengths, directly, 
Hero again, evaluation of requires the evaluation of 
One way of doing this would bo to assume the concentration 
of 13 T to be zero till tho concentration of reaches a 
(Icl’inibo vnliio. Bor examolo initially ff is taken as 10 . 

j r.j ■■ to bo zero. V’hon St ^ roaches a small value, 

^ xl ** ' 

say, 0.0001 , n is oxteiidod to 20 and 82 ^ is assumed to be 
zero and this process is repeated. The difficulty with this 
method of solution is that it takes a lot of computational 
time to solve the system of equations especially at high 
conversions. Also the truncation error or round-off error 
becomes higher, the larger being the chain length. Hence, 
recourse has been taken to tho method, employing moment 


equations. 
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i... i; 11£: considor a Goljcular weight distribution 
function ,(v), which would gonerate the kth mouiont of 

c /‘ritinuo'uo cHviin v-z3ri„\bLo vt 


\ == J (vr U{v) dY (4^1) 

0 

:i So t oi inf i.ii CG iiumbcr of nioniGnts with aii Bx^msioii 
ill wCa.*is of QT d. polynoniialG yiolds a differential 
u j, M ihn tioii * Min pointed out thst the lag^ome poly— 

i *-“1,1 u lit'ivc Qoro ganoral orthogonal properties than bho 
HcrCiito polynor.ials bocause of the possible fre.ctional order 
of distribution variable within the integrand, fhe relevant 
orthogonal properties are: 




and 


i X oxp(- x) (x) L^(x) dx = Slls^^+2l 5 ( 4 , 2 ) 

,3 .1 m nJ nm 

l“(x) is tho iu -Ihiod 'Daguerre polynomial, 

lS m =r (-i)hS!“) (4.3) 

J=b 


= 1 for n=m 

nm 

s= 0 for njka. 


(4.4) 


Any normalized distribution function &(x) can be 
reprosonted by a choice of Lagaerre expansion, as follows 

oO 


(t 

n=0 


(4.5) 
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Bivt, lor SKiArtth thu number of 

pol 5 ’'nomials x'couirod if:. :io 'c infinite. Thus the above 
rjura-'uaion io truncated at some i Ifow, using the 

orthogonal properties of Laguorre polynomials, tho coefficients 

« = J 

arc corrolatod to the moraonts of the normalized distri- 
bution r’.'.r.ctio:: G(x) by the following expression; 


as 

n 


" (-1)^ (?) a; 

rTo+o+rj 

3=0 



(4.6) 


Min [11] suggor.tod that the free parameter a bo so chosen 
that both and Cg become zero, resulting in a fast convergence 
The expression for a, according to Min [11] is 


a 


-1 + ^ 


(\)= 


( 4 . 7 ) 


Tho normalized distribution fxmetion G(x) is related 
to the original distribution function M(v) by tho following 
expresoion: 

M(v) dv = G(x)dx (4.8) 


viiLii'a represents the zoroth moment, 

(a+l) 


M- 


0 


(4.9) 


The moments of the two distributions, ond G^ are 

.ted to each other by the expression; 

~1 

s» Sf r — 2——— 


( 4 . 10 ) 
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Frov. the nbovo equations it can bo vscon that the 
dii'Xhr.'-n i; j, >1 distribution can bo .opproximatot;: fropa a finite 
nus-ibor to b.j precise) of integer mom^-nts of that distri- 

bution oy using ( 4 . 5 ) to ( 4 . 10 ). The number of such moments 
depends upon the complexity of the shape of the distribution. 

In the present work, arbitrarily chosen to be 5. Thus 

ro'ictor oodt'lling will only require ^5 before solution 

of oqua.tions can bo obtained. Moments of higher order i.e. 

6th and higher ones, can bo obtained with the help of eqns. 
(4.1) to (4.10), by the following equation: 

—It 1 

G sr f :C gCx) dx 

JL 

0 v ■ N 

max n 

«nr 

n=3 >0 

(-1)3 naty+h ,. r . (nt £!Ji^ — ^ (4.11) 

V (n-.o+i)r'(a+j+i) t^(o+i) 

To bo specific can be obtained if the first five 
mom out 8 arc known and then X g of th® original molecular weight 
distribution can be calculated from eqn.(4.lo), i.e. 


\ \ r 

>6 = Ao L-J- ^aTl) J ^ 


( 4 . 12 ) 


Further, ‘3^ ooncuntration of the amino caproic 


ncid , can be evaluated from eqn. (4.8) 


Si « 


M(v) 


Vsl 


(4.13) 
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I’ira'j D.nd g can bo exprcissod in terms of the 


Jii 


f;L, i.*-. ^ to and can bo substituted in 


o L!i’ n 

tlO' i-ysliOi oi 2 'oodolling cq’aations for the batch process. 

'fh^n .'eoduiliiic, simply constitutuc the solution of simultaneous 
'ii io rv.nti,.'.l equations and non-linear algebraic equations. 

Once the moments are calculaitcd, the numhe^r average 
v/wiglit average molecular weight, and hence the poly- 
il iaperoi ty index for the i^oljrmer product can be calculated, 

■riio O’h ■.'•■' 10 , for these quantities are as given below; 

= -v^ (4.13) 


M, 


w 




\ 

X o 

Xi 
1 


w 

f! 


(4.14) 

(4.15) 


n 


mO lu tiun pi.', HpdelJ,in£; Jqju;itipn.& f or, _tho CGTR-Pl ug Flow Reactor 
(fiubi'i.ibititi; 


'I’U..- rt..’tV4;uu of modelling eque-tions for thics combination 
ceaairliB of (;hc- norA-li’uar algebraic equations (3.13) and 
(5,15) for the CSTR, and the equations (3.3) and (5.9) for 
the plug flow reactor. 

fho per cunt by weight of water in the molt at the 
exit of the CSTR is given by the expression [lO] 


W 1.76 - 0.0060 


(3.16) 


So ofice Tq is known, the water content can be evaluated. 



27 


I’onvcr&inn oT the vr,i!;,‘r 


con'jtnt 


into liiole.r units can be 


-i' c iij. ij V o(i - tc x'xal. ^^ciancc fox' "Wcvtox 

balanco (Sc,j Appendix: l). 


.Tid the total mass 


^1 -A 6 ‘''^^ined in teitns of tbvo first 

five laomenter^s has been just indicated. iTow, eqns. (3.15) 

.,.i.nd (3.15) nro a closed set oi non— linear algebraic equations 
.■/hich can bv solved by following any suitable technique, in 
Uiu pruL’.;nfc work, the mcuhocl suggo&ted by Brown (12] hes been 
..rs.c auu of ior solving the atcove mentioned equations. The 
uloCii diagraia ior Brown’s technique is gi/en on the next pa,ge. 

rhe residence time in the CSTR was varied from 1 hour to 
5 hours, Th*-' concontratio-'S of caprolactcju, oinino caproic 
o.cid, acid end groups and the moments were computed for 
oacih value of residence time. The output stream of the CSTR 
in th- far the plug How region. Prom this point onwards, 

till th*; m'nd of the rcacbor, tho evaluation of the above 
-■'u-as li.ouvd quantities (i.o, the concentrations of the various 
especiua 'ind the momoiits) was done by integrating eqns. (3.3) 
and (5.10) by Runge Kutta method, sjid PDI were also 

c.ilculat .'d using tho moments. 



rormulo 


Gu( 2SS a solution 


olculat<2 partial 
d<2rivativ(2S 


olculotG jocobion 


/ MoximunTi partial 
V cl(2rivativ<2 


Partial pivoting based on \h<z'\ 
above to choose the variable/ 


Colculatc Lk 


Solve the equotion raving 
one unknown ond L/c 


/ Bock substitute to get 


heme tor the solution of CSTR model equotion 
(Brown's technique) 












CIIAPTijR 5 


lil - OUL'fS Aim DI3CUS3I0F 

Direct inte^'ratioa of the rate equations for (the 
concentration of chain length polymer having n units) was 
initially tried. In this method it was found that after a 
certain stage (about 5 hours of reaction time), the muaber 
of differential equations which had to be integrated were of 
the order of 600. (=n referred to on p.22 ). This consumed 
lot of time and also caused errors. For example the simi of 
the concentrations of the individual chain units was less 
than the acid end group concentr atioii (i.e.^_' <( C). 

This difference became l,..rger vrith increase in reaction time, 
rt i3 pocsiblo th.'.t duo to tiio low magnitudes of the concen- 
tratiunu, the round-off errors might have become large. Hence 
in .the present study, all the results were based on the 
nolution oi the moment oqu-ucions. 'The moments evaluated upto 
3 hours by tha direct integration method approximated well 
with those obtained by the moment equations method, thus 
proving that the second method was correct. 

Calculations were made for a CSTR followed by a PPR 
by tho method described in the previous chapter. Values of 
1,2, 5, 4 and 5 hours were selected. for the residence time in 
the CSTR, The following quantities were calculated. Concen- 
trations of the caprolaotam and acid end groups, the first 
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,.,nd thi. occond inoaonts. jpro;:i tho momonts,, the number and 


woi^,,ht. avora^o aolccul.-'.r 


ve.i.rd'its and hence the polydispersity 


indc--. .-fcro also evaluc.cod, .-ho listing of the coraputor programs 
for both tho cases have teon giv^n in Appendix II (Pago Fo.55). 
In all the ca.lcula’iions, j[ ^ been assumed that feed contains 
0.59 moles of water for every 3.75 moles of caprolactam. Those 
iif;,ures nave boon taken _prom the experimental batch data 
published by Honiians et al, [1]. This facilitates comparison 
of the model with the data of Hermans et al.fl]. 

Prom the results thus obtained, the per cent conversion 
o-i lactajn is plotted against the total residence time in the 
re. ...c tor, in i‘'ig.2. Curve j depicts this relationship for a 


(batch voac'tor). 

fec’ d ihi obi,Ain.od by Hermans ot al. [1] are also shown. 
It c lii ao .".t/.:.! cl'-.j t the rj,‘,ults predicted by the model agree 
quite w. 11 with thoO'C ODtainod by Hermans et al.fl]. Curves II 
to '/I represent tho s;rfflo relationship for a OSTR followed 
by a PPR, the variable being the residence time in the CvSTR. 
for this system the effective concentration of water in the 
CSTR (after a certain fraction of the water fed has evaporated) 
has boon determined by the method described in the previous 
chapter and is found to be 0.236 moles per kg. of the feed 
mixture to the system, it appears from Pig. 2 that the per 


cent conversion of caprolactam is maximum in the case of tbo 


batch reactor. This is so because tho water content in the 
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CSTR-PPR combination system is only 0.256 moles as compared 
to 0.59 moles in the* Pi'R. Since the concentration of water 
in the CSTR is low, the ring opening reaction (i.e. the 
hydrolysis of caprolactam) is reduced to a low level, vindi- 
cating th'-. nature of the curves in Pig. 2. fhc per cent 
conversion in 'the case of curves II to '71 exhibits an increasing 
trend with an increase in the residence time in the CSTR. 

3Jue to the strong mixing effect in the CSTR, the conversion 
of caprolactam is higher, the longer the reaction mass resides 
in the CSTR, thus explaining the nature of tLie curves II to VI 
in Pig. 2, 

In Pig. 5 the relationship between the total residence 

of 

time and the numb cr-avc, rage molocular weight ^polymer is shown. 

Prom the curves it can do seen that the numbei' average molecular 

weight of the pol^mier, in the case of bhe babch reactor 

(curve I) is greater than tho.t obtsiined in the combined system 

of reactors, in the CSTR of xfhich the residonco time is upto 

2 hours (curves II and III). But the same curves for the 

residence times of 5 hours or more (curves IV to 'Vl) show an 

improvement in the value of M over the babch system. Since 

n 

is directly proporti.onal to tho DP of tho polymer producb, 
an idea of the DP of the product can bo obtained from Pig. 5. 
Comparing the curves in Pig. 5, it is obvious that for obtaining 
the s:uflo DP, there is a rednetion in the total residence time 




Effect of tct3l residence tin^e on number overage 



34 


for the OSTR-i'I'R system with a riisidence time of 3 hours 
or more from fchat of the i-'i' i (batch system) . This can be 
explained by tho fact th:,t in the former caso, because of 
the lesser amount of wat-n"* present in tho systom, th- reverse 
condensation reaction is chccJtod to a great extent while in 
the latter (batch system) this is not so. Tho reason for an 
increase in with an increase in the residence tim.-? in tho 
CSTR is as follows: As tho residence time increases, concen- 
tration of the acid end grotips also increases, due to which 
larger chain-lengths can be obtained through the condensation ; 
reaction. If this were not so, part of the PER is utilised for 
carrying out tho poly addition reaction which does not increase : 
rapidly. 

A. similar kind of behavioui' is exhibited in the 
rol.itionsjiip bu'ewoen the overall reaction time s,nd the weight 
avorage molocnlar weight o.s shown in Pig, 4. Since weight | 

avorago uolocular weight is closely related to tho melt viscosit] 
an idea of the molt viscosity can be had from these curves. 
Jacobs and Schweigman [10] determined the effect of removal ! 

of wator at various stages of a continuous polymc^rization | 

I 

process on the relative viscosity of the polymer solution. | 

They also determined the effect on viscosity when the per 
cent removal of water from the system was varied. They conclude^ 
that the increase in the relative viscosity of the polymer i 

I 

solution was higher ^ tho Ir^tor boing tho romovo.! of water froni | 
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the system or the groatur being the per cunt romoTrnl of 
irator. 111 our system t]i>j molt viscosity of th .• polym^jr has 
boon dotorminod from instoid of M^, by adopting the method 
of Kbhan [13] . The molt viscosity was then plotted against 
the total residence time for the CSTR-PJFR system for a 

rosidonco time of 4 hours in the CSTR, as shown in Pig. 5. It 

can be seen that tho melt viscosity increases with tho total 
residence, time. This data should be of value in modelling the 
heat transfer and fluid mocha.nical aspects of tho reactor. 

A plot connecting tho total residence tine and the poly- 
disporsity index (PDI) is drawn in Pig. 6. Por the batch reactor 
the RDI gradually increases to a value of 2.0 at the end of 

the total residence tiiae. The value of 2 is ti/pical of 

poly condensation reactions. Kilkson FlAl showed that in 
a stirred tank reactoh where irreversible random condensation 
polynorization occurs, che PDI assumes very high values in the 
absence of chain terminators, while in the case of the irrevur— 
sible o,ddition polymerization reaction in a plug flow reactor, 
thu PDI reaches a finite mo.ximum value and then decreases 
gradually to a constant value. The mechanism of caprolactam 
polymerization indicates that both condensation and addition 
polymerization reactions occur simultaneously and moreov-ir , 
the reactions ai-e reversible. Also, the reactor is modelled 
to contain the CSTR and tho plug flow reactor. Thus tho PDI 
can be expected to attain a reasonably high level in the OBTR 
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sta.^:^ , and onco tho reaction mixture enters the plug flow 
region, the distribution can be* expected to got harrowed 
down to a lower level, due to the dominance of the 'reversible 
reactions, 

Tho effect of inci’casing the residence time in the CSiit 
on various quantities like the per cent conversion of lactam > 

and Pfl of the polymer is shown in Pig.?. It can be 
seen that with increase in the residence time in the CSTR, 
all the quantities increase. Out of the five values of tho 
residence time chosen in tho CSTR stage, the residence time 
of 4 hours seems to indicate the proper length of time during 
which the reaction mo.ss con be allowed to react in the CSTR. 

This conclusion is based on tho following; Tho por cent conversion 
fox’ this co,so is noarl^r tho sane as that for 5 hours residence 
tiiio in 'bho CSTR (Fig. 2), for these two cases' is also the 

sanio (fig. 3) and Phi roaches a valuo of 2.0 at the end of the 
total rosidonce time (Pig. 6). •' 

Tho rolo.tionship botwv3on the per cent Conversion of 
lactaim and (i) tho numbor-avorago molecular weight, . (ii) the 
weight-averago molecular weight and (iii) the polydispersity 
index, for the I'-eactor models considered have been shown in 
Figs. 8— 10. 'The behaviour of these curves is much similar to 
that of the curves connecting these quantities and the total 
residence time, and tho discu.ssions made for the latter curves 
hold good for the former ones also. 




conversion 
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°lo conversion of coprolactom 


Fig. 8 -Effect of °/o conversion of coprolactom on number 
rjveroqe molerulor weight 





Effect 0 
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In Pig. 10, in vliich thw por ctnt conversion of 
cnprolnctan is plotted against the PDI of the polyxier, it 
can bo sewn that curves II to VI not only approach the value, 
of 2^0 as the per cent conversion incrua,ses, but also approach 
the behaviour of the batch reactor (’curve I), for which the 
water concentration is higher. This kind of a behaviour seems 
to indicate that irrespective 'of the initial water concen- 
tration, the PDI of the pol;nner roaches the s.amc value (about 
2 . 0 ). 

The assumption of Plory’ s distribution for correlating 
the concentration of the linear dimer and li.near monomer, 
msUdo by Hoftyzer ot al. [5l and Gerdes et al. [7] was investi- 
g^/bed from the results obtained. Results are displayed in 
T<,'.blo II. It is found tiij.t during the initi.al part of the 
proc-LS, tl'n.' assumpt'ion is incorruct, while towards the end 

■Jhon CO ide.nsation rea.ction becomes predominant, the assumption 

< 

is valid. The other assumption i.e. the concentrations, of the 
linear dimer and linear monoiuv.r are equal, made by many worker 
ovei" the years, ha,s also been found incorrect duiring the 
initial stages. 
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2AB'Lii 2: COnPAHISOlT BUTyBElI m*il VALDES OL' LIIJEAR 

HOIOIiyli AllU LINEAL DlIffiR ^2 (CALCULATED 
FROM PLORY'S DISTLIBUTION) 

[Tho following tablu givue tlio values of ‘ind S 2 calculated 
fron the batch process data]. 


Time , 
hr 

Concentration of linear 
monomer, Moles/hg. mixture 

X 10^ 

Concentra.tion of 
dimer, moles/kg, 
(C-3-, r 

Sa = — ^ 

iineo.r 

,uiixti\r 

10^ 

1 

0.09256 

0.115 


2 

0,1551 

0.251 


5 

0.1729 

0.381 


4 

0.1564 

0.406 


5 

0.1563 

0.327 


6 

0.1620 

0.240 


7 

0.15 95 

0.182 


8 

0.1517 

■ 0.146 


9 

0.1417 

0.124 


10 

0.1355 

0.110 
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Gurvo I 

CurvtiS II 
to YI 


Curvo II 
Curvo III 
Curvo 17 
Curvo 7 
Curvo 71 


L-I C-GFD I-'Or. T ig. ■•'ICtIIRES 2-10 

Reprosonts tho results of tho batch reactor (PER) 
Represent the results of the CSTR followed by 
the PER with the residence time in the CSTR being 
the variable 

The variations in the residence timu in the 
CSTR are as follows; 

1 hour 

2 hours 

3 hours 

4 hours 

5 hours 



GHAP TJ3E 6 
COI CLUSIOI 

.I'wo types of reactor configuration have been considered 
in our study: the batch reactor and a stirred tank reactor 
followed by a plug flow reactor. Moment equations for the 
two systems were solved to yield the molecular weight 
distributions and the moments. It has been found that the 
PDI was in the region of ai'ound 2.0 for the s 3 ''stems considered. 

The per cent conversion of caprolactain has been found 
to be quite low in the case of the system where the stirred 
tank and plug flow reactors were combined, compared to the 
batch reactor. This has been foiind to be due to the reduction 
in the effective concentration of water in the polymer melt, 
resulting from the evaporation of water in the CSTR zone. 

The model, which appears dei'oetive due to this, can be 
corrected by assuming a closed stirred tank rea.ctor followed 
by a plug flow reactor, with water being, removed in between 
the two reactors. 

With the help of the model formulated, it is possible 
to evaluate properties of the polymer product that have hitherto 
not been calculated. An assumption for the correlation 
between the linear dimer and linear monomer that has been 
mode by various other workers, has been avoided in the 
present model. 



CHAPTER 7 


SUGGBSTIO-'^S ?OR FUTURE WORK 

The results predicted by this model could not be 
verified because of lack of sufficient data in the literature. 
Only the per cent conversion of caprolactam coulcl be compared 
with the experimental data of Hermans et al. fl]. The following 
data should be obtained for this purposes 

1. Effect of total residence time in the reactor on 
the per cent conversion of caprolactam, weighti- 
average moleculor weight, PDI and melt viscosity. 

2, Effect of varying the residence time in the CSTR 
on the above mentioned qua,ntities. 

5.V.arious sets of date, have been cited in the 
lib' 'nature for cho re-iction rate con.stants and 
!.ho correspond ill;';, equilibrium constants. Optimi- 
zatioii of those constants may be done, since the 
correct soliition of the moments depends to a 
groat extent on choosing the proper values for 
these constants. 

4. Isothermal operation assumed in this work shcild 
be relaxed. 
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APPuH SIX I 

^ 1 I TATI ■" -'AlAF: CC'M?A?T IF 
THE QgTR 


Consider the schdmatic diagram shown below. The feed 
mixture contains 8.75 moles of caprolactam and 0.59 moles of 
water. 


Hence , 


Feed 

mixture 


'I' Hater , ¥ 



put. 



I 


Mass of caprolactam in the -feed = 8 .T 52 : Mol.Wt. of 

caprolactam 
= 8.75 X 113 
= 983.8 gms 


Mass of water in the feed = 0.59 x 18 

= 10.6 g 

Total mass in the feed = ,988.8 + 10.6 = 999.4 g 
Let us assume that VI is the mass of water that 
evaporatres out of the reactor and 1 is the mass of the reaction 
mixture leaving the CSTR. 

Per cent water content in the output mixture is given 


by [10] 

H = 1.76 - 0.006 Tq 

where T^ is the temperature of the reactor 
In the present study T^ = 221°C 



Hence, 

¥ = 1.76 - 0.006 X 221 
= 0.434 pel' cent 

Drawing an overall inasB balance and a water balance 
over the CSTR, we get 
999.4 = ¥ + X 
10.6 = W + x ] 

Solving the two equations, we get 

X = 993.03 gms 


Hence, 

Water content in bhe output mixture , 
in moles 


993.03 ^ 0_,4_31 
100 Is 

= 0.236 
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